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ABSTRACT 

Context. Supemovae play an integral role in the feedback of processed material into the ISM of galaxies and are responsible for much 
of the chemical enrichment of the universe. The rate of supemovae can also reveal the star formation histories. Supernova rates are 
usually measured through the non-thermal radio continuum luminosity; however, a correlation between near-infrared [Fell] emission 
and supernova remnants has also been noted. 

Aims. We aim to find a quantitative relationship between the [Fell] at 1.26 /jm ([FeII]i 96) luminosity and supernova rate in a sample 
of 1 1 near-by starburst galaxy centers. 

Methods. We perform a pixel-pixel analysis of this correlation on SINFONI data cubes. Using Bry equivalent width and luminosity as 
the only observational inputs into the Starburst 99 model, we derive the supernova rate at each pixel and thus create maps of supernova 
rates. We then compare these morphologically and quantitatively to the [Fell] 1 26 luminosity. 

Results. We have found that a strong linear and morphological correlation exists between supernova rate and [FeII]i 26 on a pixel-pixel 
basis: 



log 



ySNn 



yr ^pc ^ 



1.01 ±0.2* log 



•41.17 ±0.9 



This relation is valid for normal star forming galaxies but breaks down for extreme ultra luminous galaxies. 

Conclusions. The supernova rates derived from the Starburst 99 model are in good agreement with the radio-derived supernova rates, 
which underlines the strength of using [Fell] emission as a tracer of supernova rate. With the strong correlation found in this sample 
of galaxies, we conclude that [Fell], 26 emission can be generally used to derive accurate supernova rates on either a pixel-pixel or 
integrated galactic basis. 
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1. Introduction 

Supernova rates (SNrates) are typically estimated by the inte- 
grated non-thermal ra dio continuum emission from their rem- 
nants (ICondonI 1 19921 and references therein). The connec- 
tion between radio continuum emission and SNrate comes 
from the tight infrared-radio relation that is usually inter- 
preted to be a natural consequence of (massive) star formation 
and stellar evolution dCondonl 1 19921 and references therein). 
Observations of near-by starburst galaxies show many com- 
pact ra dio so urces, generally attributed to supernova remnants 
(SNR). ICondo n & Yin (1990) use observations of SNR in our 
own galaxy to establish a relation between th e non-thermal ra- 
dio continuum and SNrate. iHuang et al.l d 19941) find a similar re- 
lation using observations of compact radio sources in M82. 

In addition, near-infrared (NIR) observations of SNR of- 
ten show strong [Fell] emission li ne flux at 1 257 and 1.644 
nm coincident w i th the radio peak (iGraham et al.lll987l Il990t 
lOhva et al ] 119891 Il990 t). In interstellar space, iron atoms are 
typically locked in dust grains. However, shock fronts asso- 
ciated with SNRs cause efficient grain destruction through 
thermal sputtering. This releases the iron into the gas-phase 
where it is singly ionised by the interstellar radiation field. 
In the extended post-shock region, [Fell] is excited by elec- 
tron collisions (Mouri et al. 2000), making it a strong diag- 
nostic line for tracing shocks. Since [Fell] is commonly as- 
sociated with SNRs, the [FeII]/Br'x line intensity ratio is of- 



ten used to estimat e the abundance of SNRs relative to cur- 
rent star formation. IShuII & Drain"a (il987y s shock models of 
[FeII] /Hjg imply a [FeII]/Br7 ratio of 20-70 dMoorwood & Ohval 
1988), which is in good agreement with SNR observations 
(Graham etal. 1987; Moorwood et al. 1988). Since the [Fell] 
emi ssion hfetime of a SNR is approximately 10"^ years 
(Oliva et al.ll 19891) . the [Fell] luminosity should also be a strong 
tracer of the SNrate. Qualitatively, the connection between [Fell] 
session and supernova activity is supported by a wealth of 
galac t ic and extra galactic NIR spectroscopy ( Moorwood e t al] 
119881: [Greenhouse et al. 1991; Vanzi & Rieke 1997 ) and high- 
resolution imaging (.Forbes et al.,.1993; .Greenhouse et akll 19971: 
lAlonso-Herrero et al.1120001) . which both demonstrate that [Fell] 
emission often coincides with known radio SNRs and may in- 
deed provide at least a relative measure of SNrate. 

However, the quantitative connection between NIR [Fell] 
emission and supemovae has proved to be more complex and 
a better understanding of the underlying physics is necessary. 
Radio continuum and [Fell] emission do not alwa ys co rrespond, 
even in a specific SNR. For instance, iGreenhouse et al.l (i 1997b 
found that there is little correlation between the radio continuum 
and [Fell] line emission in radio-loud SNRs in M82 and sug- 
gest that the SNR emission characteristics change as a function 
of age. In addition, Alonso-Herrero et al. (2003) obtained high- 
resolution HST [Fell] images of M82 and NGC 253 and found 
compact [Fell] emission in only 30-50% of radio all SNRs. They 
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noted that up to 73-78% of the [Fell] flux is actually diffuse 
emission, and does not originate in the discrete radio SNR which 
makes it more difficult to interpret the connection between SNR 
and [Fell] emission. Vermaas & van der Werf (in preparation) 
argue that the absence of [PII] line emission at 1.19 fim implies 
that [Fell] emission is the outcome of grain destruction by strong 
shocks. This diagnostic allows the determination of the excita- 
tion mechanism of [Fell] Vermaas & van der Werf (in prepara- 
tion) find in M83 that although the [Fell] is diffuse, the distinct 
lack of [PII] indicates that this emission nevertheless is the con- 
sequence of grain destruction by evolved SNRs which have lost 

the ir individual identiti es. 

IVan der Werf et all ([T99l . ICalzettil ([T997h and 
l^nzi &Riekel (1 1991 have attempted to correlate SNrates 
to [Fell] emission by finding a relationship between the non- 
thermal radio continuum and [Fell] flux and using a conversion 
factor to derive a SNrate. Van der Werf et al. ( 1993) uses 6 cm 
radio fluxes of NGC 6240 from Ealeset al. (1990) to estimate 
SNrates with the conversion factor from [Condon & Yin ( 1990) 
which was based on galactic SNR observations. Calzetti d 19971) 
adapted the relation fr om Van der Werf et al. ( 1993) to [FeII]i 26- 
IVanzi&Rieka ((19971) demonstrated that in blue dwarf galaxies, 
all of the [Fell] emission is plausibly accounted for by SNRs and 
presented a relationship between [Fell] luminosity and SNrates 
based on the M82 calibration. They also found a relationship 
between [Fell] and age, assuming that most of the [Fell] derives 
from SNRs. 

Thus far, there has been little agreement between the ra- 
dio continuum and [Fell] line emission in terms of determining 
SNrates. However, such a relationship is desirable because it al- 
lows us to use relatively straightforward NIR [Fell] observations 
to estimate SNrates whenever individual radio SNRs cannot be 
resolved. In this paper, we will show that a well-determined rela- 
tionship between [Fell] and SNrate exists, and explore and con- 
strain it in a sample of eleven central regions of near-by starburst 
galaxies. 

First, we will briefly introduce the galaxies in our sample in 
Section |2] and outline the observational parameters in Section [3] 
The results are discussed in Section |4] which also displays the 
spectra, Bry, Fell, and H2 line-maps and continuum maps. In 
Section |5] we explain the methods used in the analysis, and in 
Section |6] we present our analysis of the [Fell] -SNrate connec- 
tion. Section |7] explores the dependency of our results on the 
star formation history of the galaxy and Section[8]compares the 
SNrates derived from the Starburst 99 model to the more clas- 
sical radio-continuum-derived SNrates. Finally, our conclusions 
are highlighted in Section|9] 



2. Sample 

We have selected eleven bright near-by starburst galaxies at dis- 
tances ranging from 10-100 Mpc (see Table [T]). All are spiral 
galaxies with types ranging from SOa to Sd, with the exception 
of NGC 520, which is a merging galaxy. Two of the galaxies 
in our sample are seen edge-on, NGC 520 and NGC 3628. Far- 
infrared luminosities range from 1.4 X 10'% to 3.2 x lO'^L^, 
where Arp 220 with L > IQ^^Lq must be classified as an Ultra 
Luminous Infrared Galaxy (ULIRG), and NGC 6240, NGC 1614 
and NGC 7552, all with L > lO^L© are Luminous Infrared 
Galaxies (LIRGs). Our sample does not contain galaxies with 
a dominant active nucleus, but it does have four Low-Ionisation 
Nuclear Emission-line Region (LINER) galaxies and one weak 
Seyfert 2 (Sy2) nucleus. 



2.1. NGC 3628 

NGC 3628 is a near-by (~ 12.8 Mpc) HII LINER galaxy. It is 
of morphological type Sb pec and is seen edge-on. Although 
this galaxy's nucleus is dominated by a central starburst, it 
may also be h ost to a small AGN (Gonzal ez-Martm et al.ll2006l: 
iRoberts et all2001.) . NGC 3628 is a member of the Leo triplet, 
a group of 3 closely interacting galaxies. Observations of the 
HI disk reveal a very disturbed distribution, pointing to the 
interaction b etween NGC 3628 and its near-by companion 
NGC 3627 dRotsI 119781: iHavnes et al.1 Il979l) . Radio observa- 
tions show an extended radio cor e of 350 x 60 pc with a pre- 
dominately non-thermal spectrum (iHumme l 1980; C ondon et al.l 
1982: 'Reuter et al.l[T991 ). Its inner kilo-parsec is dominated by 
a significant concentration of molecular gas (iIsraeL2009.) . 



2.2. NGC 4536 

NGC 4536 is a late-type Hll/Starburst galaxy that shows [Ne 
V] line emissi on, which may indic ate that the nucleus harbours 
a weak AGN (iHughes et al.ll2005h . Previous observations indi- 
cate vigorous star formation in the central nuclear region, shown 
through Bry emission (Puxlev et al. 1988). The diffuse radio 
emission with three central peaks dVila et al] 1 19901; iLaine et al] 
2006) may indicate that this occurs in a ring around the nucleus. 
There are two ultra-luminous X-ray sources, one of which is lo- 
cated in the core o f the galaxy and possibly associated with the 
weak AGN dLiu & Bregman 2005i 

2.3. NGC 1792 

NGC 1792 is in a strongly interacting galaxy group with the 
dominant galaxy NGC 1808. NGC 1792 has surprisingly high 
luminosities, comparable to those of its more massive partne r 



NGC 1808, at many different wavelengths (Dahlem et al 



fiartner 
1991. 



Upon closer inspection, Dahlem et al. (1994) found a strongly 
asymmetric star formation distribution caused by the external 
trigger of gravitational interactions with NGC 1808. 

2.4. NGC 1084 

NGC 1084 is a HII driven Sa(s)c galaxy located ~ 16.6 Mpc 
away. More recently, this galaxy has been classified as an 
Sbc and Sb based on its B and H band images, respectively 
dEskridge et alj|2002l) . The H band image reveals a bright nu- 
cleus with an elongated bulge. Many weakly defined spiral arms 
were detected along wi th many bright knots. NGC 1084 was 
mapped at 1.49 GHz bv lCondonI (119871) . who detected a strorig 
continuum source to the south of the nucleus. Ramva et al] 
(l2007h derived a star formation rate of 2.8 M0 yr ' along with 
chaotic star formation that is not necessarily confined to the spi- 
ral arms. The northern and southern regions of the galaxy have 
different star formation rates and typical ages. The north is char- 
acterised by a series of short bursts happening intermittently over 
the past 40 Myr, possibly d ue to an interaction with a gas rich 
galaxy dRamva et al.l 120070 . The southern half of the gal axy is 
home to much younger stars, < 4 Myr (Ramva et al.ll2007h . 

2.5. NGC 1808 

NGC 1808 is the most massive member of a small group 
of galaxies. The galactic center has a super-wind, compac t 
radio sources, and a molecular ring dKotilainen et al.l Il996l) . 
IVeron-Cettv & VeronI dl985h first suspected that NGC 1808 has 
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Table 1. Summary of galaxy parameters. 



Galaxy 


RA" 


DEC " 


Morph.'' 


Activity ^ 


Dist."* 






log "Lfir 




J2000 


J2000 


J2000 




[Mpc] 


n -In 

[km s '] 


[°] 


[Lq] 


NGC 3628 


Ilh20ml7.0s 


+ 13d35m23s 


Sb pec edge-on 


HII LINER 


12.8 


847±2 


79.29 


10.14 


NGC 4536 


12h34m27.0s 


+02dllml7s 


SAB(rs)bc 


HII Sbrst 


15.4 


1802±3 


58.9 


10.17 


NGC 1792 


05h05ml4.4s 


-37d58m51s 


SA(rs)bc 


HII 


13.2 


1210±5 


62.78 


10.22 


NGC 1084 


02h45m59.9s 


-07d34m42s 


SA(s)c 


HII 


16.6 


1409±4 


46 


10.42 


NGC 1808 


05h07m42.3s 


-37d30m47s 


(R)SAB(s)a 


HII 


12.3 


1001±4 


83.87 


10.55 


NGC 520 


01h24m35.1s 


+03d47m33s 


pec 


Merger Sbrst 


30.5 


2162±4 


77.49 


10.81 


NGC 7552 


23hl6ml0.8s 


-42d35m05s 


(R')SB(s)ab 


HII LINER 


22.5 


1611±6 


23.65 


11.03 


NGC 7632 


23h22m00.9s 


-42d28m50s 


(R')SB(s)0'' 


HII 


19.3 


1535+15 


82.44 


11.43'' 


NGC 1614 


04h33m59.8s 


-08d34m44s 


SB(s)c pec 


Sbrst 


64.2 


4778+6 


41.79 


11.43 


NGC 6240 


16h52m58.9s 


+02d24m03s 


SO-a 


LINER 


108.8 


7242±45 


73.0 


11.73 


Arp 220 


15h34m57.1s 


+23d30mlls 


Sd 


LINER Sy2 


82.9 


5420±6 


57 


12.50 



Notes. 

Coordinates of galactic nucleus from NED. 

(«NED___^ 

<'^' ICorwinetan ( ll994h 

NED Hubble Flow Distance (Virgo + GA + Shapley where Hq = 73.0 ± 5 km/sec/Mpc) for NGC 520, NGC 7552, NGC 1614, NGC 6240 and 
Arp 220, metric distances for the rest. 

Heliocentric radial velocity (cz) fromradio measurement from lPrugniel et afl d 19981) except for NGC 7632 which is from optical measurements. 

Inclinatio n measured fr om Prugn iel et al.l (Il998h 

<«> Lf/R from lSandersetaL(.2003.) 

""^ Ur for NGC 7632 from lYuan et all ( 120101) 



a faint Seyfert nucleus. This has since been refuted by many 
(iForbes et al.lll99llPhillips|[T99llKrabbe et al.lll994 '). yet there 
is recent hard- X-ray ev idence to support the Seyfert nucleus 
dAwaki & Kovamalll993h . 



2.6. NGC 520 

Using long sht opti c al pho tometry and near-infrared imaging, 
[Stanford & Baicellsl ([T990l) uncovered two separate nuclei in 
NG C 520, est ablishing t his as a galaxy in the process of merg- 
ing. |Stanford& Balcelii (Il99lh went on to determine that NGC 
520 is most likely the result of a galaxy collision of a gas-rich 
and a gas-poor galaxy that happened ~ 3 x 10*^ years ago. A 
thick dust lane at PA=95° completely obscures the primary nu- 
cleus at optical wavelengths, but the secondary nucleus is vis- 
ible to the north-west of the primary. There are strong plumes 
of ionised gas, thought to repres ent a bipolar ou tflow away from 
the starburst-dominated nucleus dNorman et alJil996.) . 



2.7. NGC 7552 

NGC 7552 is a LIRG HII galaxy and host to a weak LINER. 
This galaxy has a starburst ring surrounding the nucleus. It does 
not appear to be overly disturbed, which suggests a slow bar- 
mediated star formation evolution. Due to the quiescent nature 
of the nucleus, NGC 7552 provides a near-perfect environment 
to study the bar and starburst ring around it. In addition to the 
bar, there are two dominant spiral arms an d two weaker rings 
at radii of 1.9 and 3.4 kpc (iFeinstein etaDll990l) . iForbes et all 
(Il994bl) found an inner ring with a radius of 1 kpc in the radio, 
which is not visible in the NIR continuum but is visible in NIR 
color maps and Bry emission. The dominant circumnuclear ring 
is at a distance of 850 pc from the center of the nucleus and 
the spiral arms provide a flow of molecular gas into this region 
dSchinnerer et al.lll997h . 



2.8. NGC 7632 

NGC 7632 is a member of a loose galaxy group, containing mul- 
tiple galaxies. This is a ring galaxy, which is observ ed to be dis - 
torted and bent towards the other group members (lAmlll981h . 
The interacting galaxies located south of the NGC 7632 also 
appea r strongly distorted with much absorbing material dArpI 
Il981h . 

2.9. NGC 1614 

NGC 1614 is a SB(s)c peculiar LIRG with both LINER and star- 
burst activi ty. Al though only one NIR peak has been observed, 
NefF et al] (Il990h considered the tidal tails or plumes to be ev- 
idence that NGC 1614 is the result of a merger of at least two 
galaxies. A central nucleus of 45 pc is su iTounded by a 600 pc 
diameter ring of cuiTent star formation dAlonso-Herrero et al.l 
2001). The Bry emission reveals a double peaked morphol- 
ogy (Kotilainen et al. 2001), in agreement with the radio con- 
tinuum (jC ondon et alJll982l) . Further studies of NGC 1614 by 
lOlsson et al.l (i2010!) found the molecular gas distribution dou- 
ble peaked at R = 300 pc with an additional weaker peak in 
the center of the nucleus. A radio continuum ring is also found 
at R = 300 pc which is triple pe aked, each c onsistent with the 
brightness temperature of SNR. Olsson et al.l (12010) concludes 
that the LINER spectrum is due to shocks associated with super- 
novae and therefore most likely caused by star formation and not 
AGN activity. 

2.10. NGC 6240 

NGC 6240 just falls short of being a traditional ULIRG with a 
luminosity of L/r - 10 '^L© but otherwise has all the charac- 
teristics of the class. The extended tidal tails seen in the optical 
indicate a merger. NGC 6240 hosts a double nucleus separated 
by 1 .5" + 0.1" as seen i n ground based obs ervations in the opti - 
cal dSchulz et al.lll993h and near infrared dPovon et al.lll994 . 
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yet a separation of 2" ± 0.15" as measured in [ FeII1i.64 ^rn 
(Van der Werf et al. '1993) and radio observations (ICarral et alJ 
[1990; Bales et al. 1990; Beswick et al. 2001). The high infrared 
luminosity is partially due to the pow erful nuclear starburst and 
partially to an AGN continuum (.Genzel et alj|1998h . 

2.11. Arp220 

The ULIRG Arp 220 emits 99% of its bolometric lum inosity in 
the infrared (ISoiferet al.lll984t lEmerson eta"l]|1984 . Arp 220 
is a very gas-rich galaxy, with the highest-density component 
in the nucleus, with a molecular gas density up to ~ 10 **cm"^ 
dRangwala et al.l201 lb . This high IR luminosity and intense star- 
burst activity points to a merger , supported by the ex istence of 
two nuclei, separated by ~ 1" (IScoville et al.lll997l) . Arp 220 
has an extremely hi gh extinction, with estimates running from 
A„ ~ 50 to 1 000 (ISturm et al.ll 19961: iDownes & Solomoniri998h . 
ISmith et al] (119981) resolve the radio nucleus of Arp 220 and find 
many radio point sources, which they identify as SNRs. They de- 
duce a star formation rate of 50-100 Mq yr"' and a correspond- 
ing SNrate of 1.75-3.5 yr '. 

3. Observations 

Observations of the sample galaxies were made with the 
Spectrograph for INtegral Field Observations in the Near 
Infrared (SlNFONl) at the VLT SlNFONl provides spatial and 
spectral data in the form of data cubes in J, H, and K bands. The 
SlNFONl instrument is mounted at the Cassegrain focus of the 
Unit Telescope 4 at the Very Large Telescope (VLT). 

We took observations in all three J, H, and K bands using 
a spatial pixel scale of 0.25" corresponding to a field of view 
of 8"x8" and a spectral resolution of 2000, 3000 and 4000 re- 
spectively, between October, 2006 and March, 2007. All science 
observations were taken in the ABA' -nodding mode (300s of 
object-300s of sky-300s of object), where A' is slightly offset 
from A. The object exposures are averaged during the recon- 
struction of the data cube. 

We extracted additional observations from the SlNFONl 
archive (http://archive.eso.org/wdb/wdb/eso/sinfoni/form). 
These observations include five galaxies, Arp 220, NGC 1614, 
NGC 1808, NGC 6240 and NGC 7552. As selection criteria 
we required the galaxies to be bright and near-by, and to have 
archival observations in the 0.25" spatial resolution mode in all 
three bands and 300s integration times. 

We used the standard reduction techniques of the SlNFONl 
pipeline on all observations, including corrections for flat field, 
dark current, non-linearity of pixels, distortion, and wavelength 
calibration. We obtained the flux calibration and atmospheric 
corrections from observations of a standard star. Finally, we cal- 
ibrated the continuum fluxes to match the published 2MASS ob- 
servations in the same aperture, and we determined a "flux cor- 
rection factor" for each galaxy in each observed band, and we 
applied these correction factors to the line intensities as well. 

4. Results 

4.1. Spectra 

We show the spectra of each galaxy center integrated over a high 
signal-to-noise region in the J (Figure [TJ, H (Figure |2]i and K 
bands (Figure O. The actual spatial area over which the spectra 
were integrated is over-plotted in the K band images in Figure |4] 
and|5] We have focused on non-nuclear regions to illustrate the 



spectral signature of the extended emission. The dominant emis- 
sion lines in each band are marked with dotted lines, notably 
Fell 1 26 and Pa/? in the J band. Fell i 54 in the H band and Bry, 
Hel2.o6, and multiple H2 lines in the K band. 

4.2. Continuum and line-maps 

The SlNFONl data-cubes allow us to construct J, H and K con- 
tinuum maps as well as line-maps for each of the emission lines 
detected in these bands. These emission lines provide important 
diagnostics to trace specific physical processes. The Pa/3 (1.282 
jum) and Bry (2.166 fim) HI lines trace massive, young star for- 
mation. The [Fell] emission lines, which emit most strongly at 
1.257 and 1.644 fim, are commonly used as tracers of strong 
shocks associated with supernova remnants, nuclear winds, or 
jets. There is also a wealth of ro-vibrational H2 lines through- 
out the H, and K bands, the brightest of which occur at rest 
wavelengths of 2.122 jjm (1-0 S(l)), 2.248 A^m (2-1 S(l)), and 
2.034 jum (1-0 S(2)). H2 can be excited by UV florescence from 
massive stars, or thermally by shocks from supernovae or stellar 
winds. The relative intensities of the various H2 lines indicates 
which of these physical processes is exciting the gas. 

In Figure |4] and Figure |5] we present the K band contin- 
uum, Bry, [Fell] 1 26 and H2,2.i2 line maps for each of the galax- 
ies in our sample. The continuum maps were created by find- 
ing the average continuum level at each pixel, excluding any 
emission lines. The line-maps were created using QFitsView and 
DPUSER's evaluated velocity map function "evalvelmap", de- 
veloped by the Max Planck Institute for Extraterrestrial Physics 
and available at http://www.mpe.mpg.de/~ott/QFitsView/ This 
function fits a gaussian profile to the emission line at each pixel. 

The first galaxy in Figure|4]is the edge-on NGC 3628, which 
has a strong central dust lane running from east to west. This is 
best seen by comparing the Bry and [Fell] maps. Since K-band 
Bry emission is at a longer wavelength, it is less prone to extinc- 
tion and the morphology resembles that of the K band contin- 
uum, with an additional emission peak at the western edge of the 
galaxy. The J-band [Fell] emission is at a shorter wavelength and 
is more affected by the dust lane. It shows an asymmetric hour 
glass morphology, with the brightest regions above the K band 
nucleus. The H2 emission is similar to the Bry emission, but it 
has more diffiise morphology with strong emission throughout 
the very central region. 

NGC 4536 reveals a star forming ring best seen in Bry and 
[Fell] emission. The H2 emission is concentrated at the K band 
continuum peak, with extended diffuse emission throughout the 
ring. Both Bry and [Fell] peak in the north-west corner of the 
galaxy but their peaks do not coincide. 

NGC 1792 has a K-band continuum morphology similar to 
that of NGC 4536, but the line emission is much weaker in this 
galaxy. The low signal-to-noise ratios prevent us from identify- 
ing any real structure. However, in the Bry map a small peak 
occurs at the K band continuum maximum, with structures sug- 
gesting spiral arms emanating from the nucleus. Although both 
H2 and [Fell] emission peaks at the end of these "arms", the 
actual [Fell] fluxes hardly exceed the noise level. The H2 has a 
higher signal-to-noise than the other emission lines and there are 
is distinct H2 peaks emission in the nuclear region. 

In NGC 1084, all three emission lines appear to peak on the 
K-band continuum maximum. The Bry and H2 maps show a 
somewhat similar distribution, with extended diffuse emission 
along the major axis line, and a weak secondary peak to the 
north. The [Fell] map shows a single central maximum sur- 
rounded by weak diffuse emission. 
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Fig. 1. Spatially integrated rest-frame spectra from the nucleus of each galaxy in the J band. The spatial area used in the integration 
is shown in Figure|4]and|5]by a black rectangle. Each spectrum is integrated over an 8x8 pixel area. Fell i 26 and Pay6 emission lines 
are denoted by dashed lines and marked along the top of the figure. The broad feature near 1 .25 //m in the Arp 220 spectra is an 
atmospheric artifact along with the narrow absorption features in NGC 520, and the numerous broad peaks between 1.2-1.24 and 
1.29-1.32 jt/m in the NGC 3628, NGC 1792 and NGC 1084 spectra. 
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^ig. 2. Spatially integrated rest-frame spectra from the nucleus of each galaxy in the H band. The spatial area used in the integration 
is shown in Figure |4] and |5] by a black rectangle. Each spectrum is integrated over an 8x8 pixel area. The Felli 64 emission line is 
denoted by a dashed line and marked along the top of the figure. Several spectra are degraded by residual atmospheric OH features. 
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Fig. 3. Spatially integrated rest-frame spectra from each galaxy in the K band. The spatial area used in the integration is shown in 
Figure |4] and |5] by a black rectangle. Each spectrum is integrated over an 8x8 pixel area. Bry, H2, and Hel2.o6 emission lines are 
denoted by dashed lines and marked along the top of the figure. Residual atmospheric contamination is present in the form of narrow 
peaks in the spectra of NGC 1614 and Arp 220 and as sharp absorption features in the spectra of NGC 3628 and NGC 7552. 
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NGC 1808 also has a star-formation ring, weakly seen in 
both the Bry and [Fell] maps and a few local maxima also seen 
in the H2 map. The ring appears to be very asymmetric, with 
relative bright fluxes at the eastern edge and almost no emis- 
sion from the western edge. The H2 emission is distributed more 
symmetrically around the K band peak. 

In Figure|5j we first examine NGC 520, another nearly edge- 
on galaxy that also hosts a dense dust lane running east to west. 
The Bry morphology is very asymmetric, peaking to the west of 
the K-band continuum maximum. The H2 emission peaks at the 
same location as the Bry, but has a secondary peak east of the K- 
band center. The [Fell] emission suffers from extinction by the 
dust lane and only weak [Fell] emission is seen to the south-west 
of the galaxy mid-plane. 

NGC 7552 is a face-on galaxy with a very clear star- 
formation ring around the nucleus. The Bry emission shows a 
pronounced minimum in the center , but cle ar individual max- 
ima in the ring. Unlike iForbes et al ] (ll994ah . who found the Ht 
morphology to match that of the K band continuum, we find the 
distributions of the H2 and the Bry emission to be very similar, 
with an arm or bar-like structure extending from the northern 
Bry peaks. The [Fell] emission is different: it shows roughly the 
same ring-like shape, but with emission peaks in very different 
positions. 

NGC 7632 reveals an asymmetric ring in Bry emission. The 
emission is much stronger on the east side of the galaxy in both 
Bry and [Fell] emissi on. The ea st side is in the direction of the 
other group members (l Arpll98ll) . The distortion and asymmetry 
is most likely evidence of the tidal interaction with its partner 
galaxies. 

NGC 1614 has a very compact and radially symmetric mor- 
phology in the K band, Bry and H2 emission lines. However, 
the [Fell] emission shows a crescent morphology around the nu- 
cleus. 

The double nucleus of NGC 6240 is evident in the K band 
continuum and the Bry line images. Less clearly, the double nu- 
cleus also shows up in the H2 and [Fell] images that reveal a 
more complex morphology. In all images, strongest emission co- 
incides with southern nucleus. 

Arp 220 is another compact galaxy that appears fairly sym- 
metric in the K band and all the emission line maps. The peak 
is slightly off center towards the western side of the galaxy and 
there are two separate peaks in the [Fell] emission. 

5. Methods 

To further study the correlation between [Fell] emission and 
SNrate, we have performed a pixel by pixel analysis. We eval- 
uated the continuum flux, emission line strength and equivalent 
width at each pixel in the SlNFONl data-cubes in order to gain 
insight into the variation of physical properties (such as ages and 
excitation mechanisms) across each galactic nucleus. We first 
corrected the line-maps for extinction, and used the Bry equiv- 
alent widths as an input into Starburst 99 to calculate SNrates. 
These could then be correlated with extinction-corrected [Fell] 
luminosities. 

5. 1 . Extinction correction 

In order to study the intrinsic emission line strengths, it is neces- 
sary to determine the degree to which dust obscures each nu- 
cleus. We accomplish this by identifying emission lines with 
fixed intrinsic intensity ratios, and comparing the observed to 



the intrinsic line ratio. In our NIR database, the best line pair to 
use is Pa/? (H 5-3) and Bry (H 7-4), which has a ratio of 5.88 
(Hummer & Storev 1987) over a wide range of physical condi- 
tions. Extinction is also probed by the [Fell]i.26/[Fell]i 64 ratio, 
in which both [Fell] lines ori ginate from the same upper level 
with an intrinsic ratio of 1.36 ( Nussbaumer & Storevlll988l) . For 
completeness sake, we calculated both line ratios and compared 
the extinction maps. There is little disparity in morphology be- 
tween the [Fell] derived and HI dervied extinction maps. It is 
clear, however, that the [Fell] lines are of limited use, as the 
spectral baseline they define is small and provides relatively lit- 
tle differential extinction. More importantly, the [Fell] lines have 
lower signal-to-noise ratios than the J and K lines, the [Fell] 1 64 
being particularly bad in this respect as it is degraded by resid- 
ual atmospheric OH contamination in many of our observations. 
As a result, meaningful extinction maps based on the [Fell] line 
ratio could only be derived over limited areas for a minority of 
the galaixes. Thus, the extinctions derived from the Pa/3/Bry ra- 
tio should be the more reliable and will be used to complete the 
following analysis. 

We calculated the extinction at each pixel and created an 
extinction map revealing the regions most effected by dust. 
Extinction maps for three example galaxies are shown in 
Figure |6l Table |2] lists the averaged visual extinction (Ay) over 
the observed field of view for both line ratios. The Table also in- 
cludes the integrated, extinction-corrected Bry, Pa/3, and [Fell] 
line fluxes. The maps were integrated over the full field of view 
but were filtered for low signal-to-noise pixels and pixels with 
line-widths either too narrow or too broad. The visual extinc- 
tion Ay is determin ed assum ing using the near-infrared extinc- 
tion law Aa oc A'^-^ (iMartin & Whittet.1990 ). The disparity be- 
tween the extinctions derived from the Fe and H lines mostly 
reflects the poor quality of the former, as discussed above. In the 
case of NGC 1792, the little line emission observed has a very 
low signal-to-noise ratio (cf. Figure |4]i, negatively affecting the 
accuracy of any Ay measurement, hence also the accuracy of the 
extinction corrected [Fell] and Bry fluxes. In NGC 6240, the ex- 
tremely low Bry flux makes it difficult to accurately determine 
the extinction on a pixel by pixel basis. 

In order to obtain the best possible result, we constructed 
extinction-corrected [Fell] 1 26 and Bry line maps using pixel- 
by-pixel extinctions derived from the Pa/3/Bry line ratio only. 
We also constructed pixel-by-pixel maps of the Bry equivalent 
width directly from the observed Bry; these are therefore inde- 
pendent of the assumed extinction. The Bry equivalent width and 
the Bry luminosity thus determined are the only observational 
inputs used in the calculation of SNrates from the Starburst 99 
model, described in detail in Section lS!2l below. 

5.2. Calculating SNrate 

Starburst 99 (from here on referred to as SB99) is a tool that 
mod els spectrophotometri c properties of star-forming galax- 
ies (iLeitherer et al.l 1 19991) such as spectral energy distribu- 
tions (SEDs), luminosities, equivalent widths, supernova rates 
and colours. This includes predictions of the variations in 
these properties as a function of starburst age. The mod- 
els have been calculated for 5 different metallicities (Z = 
0.04, 0.02, 0.008, 0.004, 0.001) and 3 diff'erent initial mass func- 
tions (IMF). In addition, the two extreme star-formation modes 
are considered: the continuous mode in which star formation 
proceeds continuously at a constant rate, and the instantaneous 
mode where it has the form of a delta-function starburst. 
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Table 2. Unextinction corrected Bry, Pa^S, [Fell] i 26, and [Fell] 1 64 fluxes with the Ay derived by the Pa^/Bry intrinsic line ratio. All 
measurements have a 10% error, which is the estimated calibration uncertainty. NGC 1792 has a 20% error due to the low signal to 
noise of this observation. 



Galaxy 


Pa/3 


Bry 


[Fell] 1.26 


[Fell] 1.64 






lO""* (erg s-'cm-^) 


10"''' (erg s"' cm"-) 


10"'" (erg s"' cm-2) 


10"'"(erg s"' cm"2) 


(mag) 


NGC 3628 


0.7 


1.4 


1.4 


0.6 


20.3 


NGC 4536 


12.7 


3.5 


3.4 


4.5 


3.9 


NGC 1792 


1.3 


0.3 


0.1 




2.9 


NGC 1084 


1.0 


0.3 


0.3 


0.5 


3.0 


NGC 1808 


18.3 


4.4 


9.2 


8.0 


2.7 


NGC 520 


1.2 


1.3 


1.2 


0.3 


15.3 


NGC 7552 


20.7 


6.8 


6.3 


8.1 


5.3 


NGC 7632 


1.8 


0.5 


0.5 


0.8 


3.9 


NGC 1614 


13.8 


3.5 


3.1 


3.7 


3.3 


NGC 6240 


3.3 


0.6 


8.5 


8.7 


15.6 


Arp220 


0.4 


0.5 


0.7 




15.4 



Notes. 

Av is determined assuming using the near-infrared extinction law oc /l^'-^ jMartin & Whittetll 19901) . Although the measured Ay is very high 
in some galaxies, the extinction measured in the NIR wavelength range is significantly less. 



In our analysis, we assumed a near-solar metallicity (Z - 
0.02), a Salpeter IMF (a = 235, M,o,, = lMo,M„p = IOOMq) 
and an instantaneous star formation mode. In the case of an in- 
stantaneous starburst, the model normalises the burst to an initial 
starburst mass of 10^ Mq. We calculated the average age of the 
population dominating the emission in each pixel from the ob- 
served Bry equivalent width (EW(Bry)), which is independent of 
this normalisation. SB99 provides the expected SNrate as a func- 
tion of the age thus determined. Due to the normalisation of 10'' 
IVIq, the SNrate must be appropriately scaled by comparing the 
SB99 age-dependent prediction of ionising photon flux (N(H°)) 
to the observed ionising photon flux, which in turn scales with 
Bry luminosity. This comparison yields a scaling factor that is 
directly proportional to the initial mass and initial star formation 
rate. This allowed us to determine the actual SNrates based on 
the true initial conditions of the region represented in each pixel. 

6. Analysis 

In this section, we use the observed line-maps and the SB99 re- 
sults to do a detailed analysis of each galaxy in order to better un- 
derstand the complex relationship between [Fell] emission and 
SNrate. First, we present the age and SNrate estimates for each 
galaxy. Then we will provide both a qualitative and quantitative 
comparison of the SNrates and the [Fell] luminosity. 

6.1. Age and SNrate 

Table |3]lists the Bry equivalent widths as well as the average age 
and integrated SNrate using the SB99 instantaneous burst model 
for each galaxy. Although the analysis was done on a pixel-pixel 
basis, the values given in Table[3]are either averaged (in the case 
of age and equivalent width) or integrated (in the case of the 
SNrate) over the galaxy. 

6.2. Qualitative correlation 

To illustrate the morphological relation between the [Fell] emis- 
sion and the SNrate, we show a side-by-side comparison of 
the K-band continuum, the extinction-coiTected Bry flux, the 
extinction-corrected [Fell] 126 flux, and the derived SNrate in 
Figurelllfor the galaxies NGC 4536, NGC 1808, and NGC 7552, 



Table 3. Average equivalent width of Bry and the SB99 derived 
average ages and integrated SNrates using the instantaneous star- 
burst model. 



Galaxy 


EW(Br7) 


Age/„„ 


SNrate,„„ 




A 


(Myr) 


yr"' 


NGC 7632 


5.9 


7.0 


0.03 


NGC 3628 


7.5 


6.9 


0.01 


NGC 4536 


10.6 


6.7 


0.3 


NGC 1792 


4.1 


7.9 


0.007 


NGC 1084 


4.8 


7.3 


0.009 


NGC 1808 


7.2 


6.9 


0.06 


NGC 520 


11.7 


6.6 


0.3 


NGC 7552 


13.7 


6.7 


0.3 


NGC 7632 


5.9 


7.0 


0.03 


NGC 1614 


22.9 


6.4 


0.9 


NGC 6240 


3.6 


7.7 


3.6 


Arp 220 


8.0 


6.8 


0.7 



which have the highest-quality observations. A visual compari- 
son of the morphologies in the four different maps of each galaxy 
shows that the SNrate map most closely resembles that of the 
[Fell] emission. 

For instance, in NGC 4536 there is a bright knot of Bry emis- 
sion directly north of the nucleus with secondary emission peaks 
to the north-west and south-east of the galaxy. The [Fell] emis- 
sion is concentrated on the galaxy nucleus with a long plume of 
emission (perhaps an inner spiral arm) extending from the nu- 
cleus towards the north-west. A small knot of faint [Fell] emis- 
sion coincides with the Bry peak. Mimicking the [Fell] emis- 
sion, the SNrate peaks at the center and also shows an arm ex- 
tending towards the north-west. The SNrate exhibits a slight in- 
crease at the Bry peak, but on the whole more closely resembles 
the morphology of [Fell] emission. 

The starburst ring surrounding the bright nucleus of NGC 
1 808 shows up clearly in Bry. The ring contains individual bright 
knots of Bry emission. The ring is much less prominent in [Fell] 
emission. The SNrate map shows a ring structure more clear than 
in the [Fell] map but it is quite diffuse and lacks the contrast that 
the Bry knots provide. The SNrate map also shows diffuse ex- 
tended structure around the nucleus, similar to the [Fell] emis- 
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sion but distinct from the Bry map in which the compact nucleus 
is much more isolated. 

Finally, the face-on spiral, NGC 7552, also reveals its star- 
burst ring in Bry, with little emission at the nucleus. In the north- 
ern part of the ring we find three major Bry knots, two slightly 
north-west of the nucleus and one elongated knot directly north 
of the nucleus. Comparison shows that only one of the two 
north-western Bry knots is discernible in [Fell] emission. We 
also see that the morphology in the elongated knot differs in 
the two maps. Focusing on the SNrate morphology, the north- 
ern emission knots resemble the relative flux and morphologies 
of the [Fell] emission much more closely than the Bry emis- 
sion. Specifically, the top north-western knot is also missing in 
SNrate map and the elongated northern Bry peak is resolved to 
two peaks, resembling the relative flux ratios seen in [Fell]. 

All comparisons show that, at least qualitatively, the SNrate 
maps resemble the [Fell] maps much more than those of the Bry 
emission on a pixel-by-pixel basis. Once again we emphasise 
that the only observational input into the SNrate calculation is 
the Bry luminosity and EW(Bry); the [Fell] flux is never used. 
As the SNrate nevertheless correlates better with the [Fell] than 
with the Bry emission, we conclude that this strongly supports 
[Fell] as a robust tracer of SNrate. 



6.3. Quantitative correlation 

Now that we have established that a qualitative correlation exists 
between [Fell] emission and SNrates, it is important to verify to 
what extent this relationship holds up quantitatively. To do this, 
we compare in Figure|8]for each galaxy the [Fell] 1.26 luminosity 
and the SNrate derived from Bry pixel-by-pixel. The line-maps 
were not additionally filtered, with the exception of NGC 6240 
which was filtered to eliminate pixels where the Bry flux is insuf- 
ficient to determine the true extinction. We normalised both the 
SNrate and the [Fell] luminosity to values per square parsec so 
that each pixel represents the same physical area. We performed 
a least-squares regression in log space on the combined points 
from all sample galaxies except NGC 1792 because of its poor 
signal-to-noise data. The combined linear regression is indicated 
by a solid black line in Figure[8j it takes the analytical form: 

log^— — -z - (1.01 + 0.2) * log I , -41.17 + 0.9 (1) 



36 



yr ^pc ^ 



ergs ^pc ^ 



The errors represent the standard deviation of the slopes and 
intercepts of the individual galaxies, thus representing the vari- 
ation between individual galaxies' fits. This is an intrinsic error, 
corresponding to the variations in properties among the galaxies 
and not derived from the quality of the observation. To investi- 
gate the strength of the relationship on a global basis, we have 
plotted , excluding NGC 1792, the integrated SNrate plotted 
against the integrated [Fell] luminosity for each galaxy, again 
with the best-fit regression line. Errors in the extinction are a 
potentially significant source of uncertainty in the derived rela- 
tionship. The dashed lines at either side of this best fit in Figure|9] 
mark the regression lines we would find if the extinction magni- 
tudes were overestimated respectively underestimated by factors 
of two. It is obvious from the Figure that our result is essentially 
unaflfected by this uncertainty. The best fit for the global values 
is given by: 
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Fig. 8. A pixel-pixel plot of SNrate, as derived from SB99, com- 
pared to [Fell] luminosity. Each galaxy is represented with a dif- 
ferent color and the values are normalised to a square parsec. The 
black line represents the best fit power law excluding NGC 1792, 
which is represented by diamonds. 



Thus, the averaged and the integrated SNrate - [Fell] lumi- 
nosity relation are identical within the errors. Both procedures 
yield a power law with a slope of nearly unity. The relation is 
linear well within the errors. 

An important source of uncertainty in the SNrate calculation 
is caused by observational errors in the Bry, [Fell] and K band 
continuum intensities. These observational errors are dominated 
by the uncertainty in the calibration, which is of the order of 10% 
of the flux in all cases except NGC 1792, where the errors are 
dominated by noise, estimated to be about 20% of the flux. The 
uncertainties in the derived SNrate are the same as those in Bry 
flux, i.e. about 10%. However, the observational errors are not 
the dominant source of uncertainty. This is the systematic error 
introduced by the choice of the starburst model. Because we use 
SB99, our choices are limited to two star-formation models only, 
instantaneous and continuous. In the next section, we discuss the 
effects of the choice of the burst model on the relation between 
SNrate and [Fell] luminosity. 

7. Dependency on burst model 

So far, we have based our analysis on the assumption that star 
formation in the sample galaxies is best represented by an instan- 
taneous burst. The other extreme case also modelled by SB99 is 
the continuous star-formation scenario. It is unlikely that either 
of these extreme cases is a true representation of the situation in 
the galaxies considered; we suspect that reality is somewhere in 
between and better described by one or more star bursts extended 
in time. We cannot determine the actual star formation history of 
the sample, but the SB99 model allows us to constrain the valid- 
ity of our determination of SNrate/[Fell] luminosity relation by 
performing the same analysis this time assuming a continuous 
star-formation mode. The result is shown in (FigurefTOb as a plot 
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Fig. 9. Integrated SNrate plotted against integrated [Fell] lumi- 
nosity over the field of view. The black line is the least-squares 
regression line, excluding NGC 1792. The red dotted lines repre- 
sent the least squares regression line if the magnitude of extinc- 
tion was overestimated respectively underestimated by a factor 
of two in optical depth. 

of the [Fell] luminosity versus the SNrate for the combined pix- 
els of all sample galaxies. In this Figure, the black dots represent 
the SNrates assuming an instantaneous burst (as in Figure[8]) and 
the red points represent the same pixels where the SNrate is cal- 
culated assuming continuous star formation. 

From Figure [TO] it is clear that also in continuous star- 
formation mode, the SNrate is still closely related to the [Fell] 
luminosity. However, in this case, the relation is no longer linear 
as the best fit slope is 1.20. We also note that in this situation, a 
given [Fell] luminosity also implies a SNrate roughly an order 
of magnitude higher than in the instantaneous burst model. 

8. Comparison to Radio SNrate 

Which mode best represents our sample? Radio continuum emis- 
sion is the classical tracer of SNrate, and radio flux densities 
can be used to obtain an independent estimate of the SNrate for 
comparison to the SNrates derived with t he SB 99 models. We 
apply the equation given by iHuang et al.l (Il994t) . which relates 
non-thermal radio luminosity and SNrate, to VLA observations. 
There are no VLA archived observation of NGC 7632 and this 
galaxy is excluded from the comparison. A similar situation ap- 
plies to NGC 7552, but here we could use flie ATCA 4.8 GHz 
map instead. For the nearest extended galaxies, we matched the 
radio integration area to that of SINFONI. For the more distant 
NGC 1614, NGC 6240, and Arp 220, we used NED VLA 1.4 
GHz integrated flux densities since these galaxies are compact 
enough that the entire galaxy is encompassed in the SINFONI 
field of view. In Figure [TT] we compare the radio SNrates thus 
derived to the SB99 instantaneous and continuous SNrates. 

For all galaxies, except NGC 1792, NGC 4536, NGC 6240 
and Arp 220, the radio SNrate closely matches those given by 



Fig. 10. Pixel-pixel comparison of SNrate to [Fell] luminosity. 
The black points represent the pixel values for all the galaxies 
using the instantaneous-burst model, while the red points repre- 
sent the same pixels using the continuous-burst model. 
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Fig. 11. Comparison of the SB99 instantaneous SNrate (x-axis) 
to the radio SNrate (blue diamond), the SB99 continuous SNrate 
(red asterisk) and for comparison the SB99 instantaneous SNrate 
(black asterisk). 



the SB99 instantaneous-burst model. As noted before, the NIR- 
emission line-maps of NGC 1792 are too noisy to be reliable. 
However, NGC 4536 is a nearly face-on galaxy with relatively 
low extinction and high signal-to-noise spectra, and the NIR data 
should be reliable. Thus, NGC 4536 is undergoing star formation 
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that is closer to instantaneous but extended in time. Arp 220 and 
NGC 6240 have radio SNrates even higher than the SB99 instan- 
taneous SNrates. Also, in Figure [8] the points for these galaxies 
lie systematically above the best fit line, with a similar slope to 
the rest of the galaxies. This offset can be explained by either an 
excess of [Fell] emission or a deficiency of Bry emission, since 
these are the two observables that are used to scale SNrate and 
[Fell] luminosity. In the case of Bry deficiency, the Bry emission 
could be selectively suppressed by dust absorption of Lyman 
continuum photons in these very dense systems. However, in 
NGC 6240 and Arp 220, it is tempting to speculate that this ex- 
cess is caused by merger-related shocks that are over-exciting the 
[Fell] in addition to the SNRs. Evidence for large-scale shocks 
in NGC 6240 comes from the prominent w ings on the H2 1-0 
S(l) line which trace the shocked superwind dVan der Werf et all 
1199 3: Engel et al. 2010). The [Fell] emission is morphologically 
similar to the H2, indicating that the [Fell] also traces the su- 
perwind. These powerful merger (U)LlRGs may demonstrate a 
limit to the validity of the direct relationship between SNrate 
and [Fell] luminosity. In addition, determining an accurate ex- 
tinction in (U)LIRGs presents a challenge. The under or over 
estimation of extinction leads to inaccurate [Fell] luminosities 
as well as Bry luminosities used to scale the SNrate. However, 
the offset between the (U)LIRGs and the other galaxies is mi- 
nor They share a common slope and appear to be offset only 
in [Fell] emission. The validity of this relationship in the case of 
(U)LIRGs can only be established by studying a larger sample of 
these galaxies. In any case, it appears that the majority of (mod- 
est) starburst galaxies is well-represented by the assumption of 
a (nearly) instantaneous burst of star formation. In addition, the 
very good agreement between SNrates derived from the radio 
continuum, and from NIR data provides added confidence in the 
diagnostic strength of [Fell] as a tracer of SNrate. 

9. Conclusion 

Using SINFONI observations of 11 near-by galaxies, we have 
performed a pixel-by-pixel analysis of the correlation between 
[Fell] 1 26 emission and independently derived SNrates. We mea- 
sured accurate Bry, Pa/3, [Fell] 1.26, [Fell] 1 64, and H2,2.i2 hne 
fluxes. We determined Bry equivalent widths which were used 
as input into the SB99 model to find the starburst age in each 
pixel. In the SB99 model, this defines a SNrate, normalised to 
an initial mass of IO^Mq, which we scaled with the observed 
Bry luminosity to find the actual SNrate in each pixel. 

The comparison of the [Fell] 126 luminosity to the SNrate 
(derived only from Bry equivalent width and luminosity) reveals 
a nearly linear correlation, where the error represents the spread 
of intrinsic properties over our sample of galaxies. 

l^g^^SNrare ^ ^ ^ ^ ^ ^ to^i^^^^^^ - 41.17 ± 0.9 



yr ^ pc ^ 



ergs ^ pc ^ 



This relationship is valid both on a pixel-pixel basis and for 
the integrated galaxy. For the integrated [Fell] luminosity and 
SNrate, the fit is remarkably tight with very little spread. 
However, to correctly determine the absolute SNrates, it is still 
critical to determine whether star formation has occurred in 
a (nearly) instantaneous burst, or has proceeded in a continu- 
ous fashion. SNrates derived from radio continuum observations 
may be used to distinguish these scenarios. Most of the modest 
starburst galaxies in our sample are best fitted assuming instan- 
taneous star formation. However, we find that the relationship 
breaks down for (U)LIRGs in our sample. 



From the strong correlations found in this sample of galaxies, 
we confidently conclude that [Fell] 126 emission line strengths 
are quantitatively correlated with SNrates and can be used to de- 
rive these rates either locally or globally. This use of [Fell] as a 
robust tracer of SNrate provides us with a very useful diagnos- 
tic tool. It is particularly important as it allows determination of 
SNrates from NIR observations of distant galaxies, where indi- 
vidual SNRs can not be resolved. 
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K Bond Continnum Bry [Fell], 26 




Fig. 4. K band continuum map, Bry, [Fell] 1.26, and Hj.i.n line maps of NGC 7632, NGC 3628, NGC 4536, NGC 1792, and NGC 
1084. The pixel scale is given by the 2" scale bar in the K band column. The black rectangle marks the areas over which the spectra 
in Figures [ri|2]and|3]were integrated. In all figures, north is up. 
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K Band Continnum Bry [Fell], je H2 2 




Fig.5. K band continuum map, Bry, [FeII]i.26, and H2,2.i2 line maps of NGC 1808, NGC 520, NGC 7552, NGC 1614, NGC eil^ 
and Arp 220. The pixel scale is given by the 2" scale bar in the K band column. The black rectangle represents the area over which 
the spectra in Figures [T]|2]and|3]were integrated. In all figures, north is up. 



Rosenberg, M. J. F. et al.: [Fell] as a tracer of Supernova Rate 





A, (mog) 



Fig. 6. Extinction maps of NGC 4536, NGC 1808 and NGC 7552. The color scales are displayed to the right of each map, the values 
represent Ay measured in magnitudes. 
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Fig. 7. A comparison between the K-band continuum, Bry line map, [Fell] i 26 line map, and the SNrate morphologies for NGC 
4536, NGC 1808, and NGC 7552. The black discontinuities represent areas that have been filtered due to low signal-to-noise values 
in the line or extinction maps. The white pixels represent the highest flux regions. 
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